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The underlaying assumptions and uncertainties involved in the derivation of WIMP exclusion 
limits from XENON10 and XENON100 detectors are examined. In view of these, recent claims 
of sensitivity to light-mass Weakly Interacting Massive Particles (WIMPs) are shown to be over- 
stated. Specifically, bounds constraining regions of interest in WIMP parameter space from the 
DAMA/LIBRA, CoGeNT and CRESST experiments can be assigned only a very limited meaning. 

PACS numbers: 



The possibility of relatively-light (<10 GeV/c 2 ) 
Weakly Interacting Massive particles (WIMPs) being the 
source of the annual modulation effect [l[ observed in the 
DAMA/LIBRA experiment 0] was first proposed in Q 
and later revisited in [4| . Since then, it has gained in pop- 
ularity as constraints from other searches have depleted 
alternative dark matter scenarios, and several particle 
phenomenologies have generated plausible candidates in 
this mass range. The situation has recently gained in 
complexity with the observations from CoGeNT [H, [f| 
and CRESST Q, which may point at a light- WIMP pa- 
rameter space compatible with DAMA/LIBRA. 

Studies of detector sensitivity to such light dark 
matter particles must be regarded with caution Q. 
For presently-available technologies, light- WIMP signals 
would fall uncomfortably close to detector thresholds, a 
region where systematic effects can lead to rushed claims 
of exclusion or detection. Experimentalists should at- 
tempt not to aggravate what is a naturally murky area 
of study, by describing the assumptions made to arrive 
to their conclusions, and by including an assessment of 
the uncertainties involved in their analyses. Recent work 
by the XENON10 and XENON100 0, EI collabora- 
tions is lacking in these respects. The goal of this report 
is to provide this missing information as succinctly as 
possible. 



I: RECENT MEASUREMENT OF £ eff BY G. 
PL ANTE ET AL. [Tlj. 

A previous attempt to derive light- WIMP limits by 
the XENON100 collaboration [12| was received with crit- 
icism [l3l | pointing out the unphysical behavior of the 
quenching factor in the production of direct scintillation 
by nuclear recoils (£ e ff) employed. An additional cri- 
tique Q emphasized traceable mistakes made in the data 
analysis of previous £ e ff measurements performed by the 
XENON10 and XENON100 collaborations. 

A new recent £ c s measurement by the XENON100 col- 
laboration (Plante et al. [TH) addresses the concerns in 
[H, specifically the systematic effect introduced by erro- 



neously normalizing simulated recoil rates to their cor- 
responding measured values, and the sub-optimal design 
of several earlier calibration detectors, prone to multiple 
scattering involving inert materials. The new detector 
used in ll| features a compact design that bypasses this 
concern and maximizes light collection from the active 
liquid xenon (LXe) volume. Not surprisingly, the mono- 
tonically decreasing behavior of £ e ff towards zero energy 
predicted in [1, [l3| is now observed by Plante et al. 

While great strides towards a better understanding of 
£ fi have been made in ll| , significant room for improve- 
ment remains: 

• An unnecessary degree of freedom in the fits comparing 
LXe scintillation measurements and simulations has been 
introduced by Plante et at, namely the energy resolution 
as a function of recoil energy, which is a predictable quan- 
tity, and not independent of £ e ff, as implicitly postulated 
in . This is in contrast to an earlier measurement by 
Manzur et al. 14| , also correctly pointing at a decreasing 
£ c g, where the resolution was determined by measure- 
ments at energies well-above any threshold effects, and 
for all lower energies unambiguously defined according to 
its expected dependence on photoelectron yield (a func- 
tion of £ c ff ) . The introduction of this gratuitous degree 
of freedom can reinstate the deleterious effects described 
in [H, by substantially biasing C c s towards artificially 
large values and reducing uncertainty. This concern is 
particularly important below ~6.5 keV r , where thresh- 
old effects become dominant in ll| . The approach taken 



llj does not necessarily have to constitute an issue, 



as long as the obtained best-fit resolution follows the ex- 
pected behavior 1 . No mention of this comparison is made 
in As described in Sec. II below, the extrapolated 
behavior of £ c g to zero energy critically determines LXe 
sensitivity to light- WIMPs, making attention to such de- 
tails very important. A discussion of this comparison 
between expected and best-fit energy resolution would 



1 Unfortunately, a sizable mismatch seems to be involved [IB 



2 



considerably improve the credibility of the lowest-energy 
£ e ff values obtained by Plante et al.. 

• Measurements in [ll[ were performed in single-phase 
mode, i.e., in the absence of the electric drift field present 
during the operation of the XENON100 detector. This 
field is expected to suppress electron-ion recombination, 
reducing the scintillation yield. While this effect was 
found to be small by Manzur et al. [3] , the C c g values by 
Plante et al. should be considered an upper limit to the 
actual nuclear recoil scintillation yield in the XENON100 
detector. This consideration as an upper limit is revisited 
in Sec. II within a different (instrumental) context. 

• It must be kept in mind that the definition of C c g used by 
the LXe detector community differs from the traditional 
one for a quenching factor, by relativizing the scintillation 
yield from low-energy nuclear recoils to that from elec- 
tron recoils at a relatively high ionization energy (122 
keV). The more conventional definition uses the ratio 
of scintillation yield from nuclear and electron recoils of 
identical energy. This may seem like a moot point, until 
the large non-proportionality ty pica lly observed in heavy 
scintillators [16| . including LXe [171 ] ■ is examined: a large 
increase in scintillation yield for electron recoils (the de- 
nominator in the traditional definition of quenching fac- 
tor) is typically observed below few hundred of keV down 
to few keV. Compton scattering measurements comple- 
mentary to those in ll| are clearly advisable. 
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II: 



NEW LIGHT- WIMP LIMITS FROM 
XENON100 [10] 



The analysis of a 100 day exposure from the 
XENON100 detector [l(| has resulted in a claim of size- 
able improvement in light- WIMP sensitivity with respect 
to a previous shorter (11 day) run 12j. A discussion of 
the strong assumptions implicitly made to arrive to this 
conclusion and of the uncertainties neglected in [Io| is 
provided below. 

• Contour "a" in Fig. 1 is similar to the exclusion curve 
in [lo| . It is obtained by assuming the logarithmic ex- 
trapolation of the £ e g from Plante et at, as proposed in 
[lOj . and that only the lowest in energy of the three ac- 
cepted nuclear-recoil events in could be due to a light- 
WIMP. Contour "b" in the same curve represents the ex- 
clusion obtained when the 2a C.L. uncertainty band in 
this C c s is adopted instead. The resulting change in sen- 
sitivity is much larger than what is indicated by the very 
narrow uncertainty bands in Fig. 5 of [1(|. This issue can 
be confirmed by performing a self-consist ency test be- 
tween the XENON100 exclusion curves in [Hf and [13: 
the two values of C c s contemplated in [H] (Fig. 1 there) 
generated exclusion curves diverging by a very large fac- 
tor for light WIMP masses (Fig. 5 in 12j). Those two 
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FIG. 1: 90% C.L. XENON100 exclusion contours obtained 
under the assumptions discussed in the text. Dotted (black) 
curves correspond to the £ e s by Plante et al. [ll[, dashed 
(blue) to that by Manzur et al. [14| , and dash-dotted (red) to 
a most recent Monte Carlo-independent C c s determination by 
ZEPLIN-III 37j . The notation used to describe each case lists 
number of irreducible recoil events accepted, £ e ff considered 
(central value of logarithmic extrapolation or its lower C.L. 
boundary), and statistics used (see text): a) 1 event, central, 
Poisson; b) 1 ev., 2a, Poisson; c) 4 evs., 2a, Poisson; d) 4 
evs., 2a, binomial; e) 1 ev., central, Poisson; f) 1 ev., central, 
binomial; g) 1 ev., la, Poisson; h) 4 evs., la, binomial; i) 1 ev., 
central, Poisson; j) 1 ev., central, binomial; k) 4 evs., central, 
binomial; 1) 4 evs., la, binomial. Additional instrumental 
uncertainties not reflected in this figure are listed in the text. 



values of C e s are coincidentally not very different from 
the central and 2a C.L. boundaries of the C e g from Plante 
et al. (Fig. 1 in [10]). However, in Fig. 5 of 0, the new 
XENON100 analysis assigns an insignificant impact on 
the exclusions to this large spread in C g. The origin for 
this lack of self-consistency must be addressed 2 . 

• In a departure from the blind analysis initially intended 
by the XENON100 collaboration, three events next to 
threshold were rejected immediately following unblinding 
10, 19]. These events have been ascribed to photomul- 
tiplier (PMT) noise affecting the SI (direct scintillation) 
channel. Post-unblinding corrective actions are often re- 
quired and no judgment on this decision should be passed 
until more details become available. However, it is worth 
remembering that this type of PMT noise was already 
present in XENON10 data [l8| (and not rejected a pos- 
teriori) and is also ubiquitous in a XENON100 example 
event catalogued as "good" [20( , indicating that data cuts 
originally deemed as adequate must have been in place 
against it 3 . Details about the post-unblinding criteria 



This question may be extended to higher WIMP masses. For 
those, following basic statistical estimators, only a marginal in- 
crease in 90% C.L. sensitivity should be expected in going from 
zero to three irreducible events following an increase in exposure 
by a factor of ten. A much larger gain in heavy WIMP sensitivity 
has been claimed in going from [12| | to [ljj. 

A single example of noise-correlation provided in [2(ill corresponds 
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developed to reject these events while in the presence 
of robust S2 (ionization) signals will be of special inter- 
est. Contours labelled "c" , "d" and "h" in Fig. 1 display 
the non-negligible effect of including these three rejected 
events into the calculation of XENON100 exclusions. 

• Dashed lines in Fig. 1 represent XENON100 exclusions 
using a logarithmic extrapolation of the alternative £ e ff 
obtained by Manzur et al. [14j using another optimally- 
designed LXe chamber. Differences in data treatment 
and mode of operation between [llj]and [3 are described 
above. Criticisms concerning the data treatment in [l4| 
are given in Q. Contours labelled "g" and "h" use the 
logarithmic extrapolation of the lower la C.L. boundary 
rather than the central £ e ff value (contours "e" and "f " ) . 

• A lingering critical question is to what extent a deter- 
mination of £ e H performed using highly-optimized com- 
pact calibration detectors like those in [ll|, Hi can be 
applied with confidence to a much larger device like the 
XENONf 00 detector ^featuring a small SI light-detection 
efficiency (just ~6% |2l|). different hardware trigger con- 
figuration, data processing, etc. For instance, simulations 
like those used within XENON100 to obtain a cumula- 
tive cut acceptance near threshold can only be regarded 
as best-effort estimates. Their limited meaning and ten- 
dency to significantly overestimate near-threshold effi- 
ciency has been recently encountered by Plante et al. 

ll| , even for the near-ideal conditions of their small cal- 
ibration chamber (~18 c.c. of LXe, with Att PMT cov- 
erage). Another example of instrumental constraints is 
the negligibly small low-energy effective £ e ff derived by 
the f 2 kg ZEPLIN LXe dark matter detector [H] when 
its measurement is attempted in situ (if adopted, this 
£ e ff generates essentially no LXe sensitivity to WIMP 
masses below 10 GeV/c 2 [Hj]). Going back some time, a 
dramatic deficit in observed neutron-induced recoil rates 
compared to few-keV r expectations was observed with 
the XENON10 detector [23(. Such comparisons should be 
revisited within the context of the existing XENON100 
neutron calibrations: if the expected response to few- 
keV r recoils is still absent due to instrumental limita- 
tions, light- WIMP limits should not be distilled from 
sheer wishful thinking. As discussed next, light- WIMP 
limits are obtained by XENON100 under the strong as- 
sumption of positive Poissonian fluctuations in scintil- 
lation light from WIMP-induced recoils well-below the 
detector energy threshold. The agreement between ex- 
pected and observed neutron-induced recoil rates should 
therefore be demonstrated into that energy range. The 
XENON100 collaboration is invited to produce this much 



needed validation of their claims 4 . 

• The mentioned assumption of Poissonian statistics gov- 
erning the microscopic processes 5 of light and charge gen- 
eration by few keV r nuclear recoils in LXe is not only 
presently unwarranted [3] , but seemingly counter to the 
scarce available experimental information: the very small 
value of the Fano factor in LXe is for instance indica- 
tive of sub-Poissonian statistics ruling those processes 
26]. Similarly, the electron emission statistics by few- 
keV, heavy-mass ions during surface collisions is known 
to be better described by binomial rather than Poisso- 
nian statistics 27]. Contours "d", "f" and "h" in Fig. 1 
display the effect of a small deviation from the Poisson as- 
sumption (a binomial distribution of same mean, taking 



a probability of SI photon detection of 6% |2l|). These 



contours should be considered as illustrative ansatzes for 
information-carrier statistics that could generate even 
less light production. In this respect, it is worth empha- 
sizing that a 7 GeV/c 2 light- WIMP is expected to impart 
a mean recoil energy in LXe of just ^0.6 keV r , with an 
absolute maximum (occurring with infinitesimal proba- 
bility) of ~4.1 keV r . At this very endpoint, the prob- 
ability of surpassing the XENON100 four-photoelectron 
(~8.5 keV, ) threshold is never larger than ~10%, even 
when the logarithmic extrapolation of the £ e ff by Plante 
et al. and Poisson fluctuations are adopted. Not all light- 
WIMP detecting media fare as poorly from this point of 
view of generation of information carriers: the same 7 
GeV/c 2 WIMP at its spectral endpoint (~L4 keV ion- 
ization) in CoGeNT germanium diodes would generate a 
readily detectable ~470 electron-hole pairs. 

Fig. 1 includes the present uncertainty in the quench- 
ing factor for sodium recoils in DAM A/LIBRA (l3[|28Tj. 
a subject of discussion avoided by the XENON100 col- 
laboration in [ijj H3|. This uncertainty extends the 
DAM A/LIBRA region to considerably lower WIMP 
masses than what is represented in 

ESQ. 



Ill: LIGHT- WIMP LIMITS FROM XENON10 VIA 
IONIZATION SIGNALS d 



A recent reanalysis of XENON10 data uses strictly the 
ionization channel in that detector to impose limits on 



to an event well-below threshold (3.1 SI photoelectrons) and does 
not seem to display the SI region corresponding to the S2 pulse. 



4 This request for validation is not unfair: COUPP results consis- 
tently include a sensitivity penalty whenever any disagreement 
with neutron-induced recoil rate expectations arises I24|l . A Co- 
GeNT detector presently installed at Soudan has not been ex- 
posed to neutron sources to avoid activation during a search for 
an annual modulation [f|. However, the response to sub-keV r 
nuclear recoils was measured with a detector identical in prop- 
erties (energy resolution, noise, bias) and crystal mass l'2. r | . 

5 Notice reference is made to those ab initio processes and not the 
ensuing statistics of photoelectron generation in PMTs. 
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light-mass WIMPs [§]. This S2 light emitted via electro- 
luminescence from charge drifted into the gas phase of the 
device is, when examined alone, sufficient to extend the 
sensitivity of LXe detectors to recoils of 0(1) keV r . Not 
including the information from SI (direct scintillation) 
allows a reduction in threshold at the expense of losing 
the ability to distinguish nuclear from electron recoils. 

While this approach is promising, the few keV r nuclear 
recoil energy scale corresponding to this S2 channel is 
presently hopelessly ill-defined. This is a result of the 
inadequate "best-fit Monte Carlo" method [29[ employed 
to arrive at it. An extensive critique of this method, 
ignored thus far by XENON10 authors, can be found in 
[Sj. In a troubling case of double-standards, the gist of 
this critique (that with this method all uncertainties are 
absorbed into the energy scale) has been recently echoed 
by the XENON100 collaboration (Sec. I in [Tl|), when 
rebuking indirect measurements of £ c g using the same 
methodology. 

In lieu of a reiteration of the criticisms in 18], the 
reader is invited to inspect Fig. 2: each of the colored en- 
ergy scales shown there, generated by the "best-fit Monte 
Carlo" method, has been claimed to be the correct one 
by the XENON10 collaboration over the brief span of the 
last two years. The scale is observed to change as rapidly 
as from workshop presentation to its published proceed- 
ings. Its monotonic evolution has been towards the black 
curves, held in Q to correspond to the most plausible en- 
ergy scale (one derived from an earlier method laid out 
by XENON 10 authors: see pertinent discussion in [|| ) . A 
critically-minded reader would (rightly) argue that none 
of these can be presently assigned any credibility at few 
keV r . However, reasons have been provided in |8j to 
support the solid black curve, representing the Lindhard 
theory modified below a ^40 keV r kinematic threshold 
[13] by an example of adiabatic correction, as in [3(3]: 

• This energy scale generates a similar quenching for the 
ionization yield and the £ e ff observed by Plante et al. or 
Manzur et al., i.e., a monotonic decrease in the genera- 
tion of information carriers (free charge, direct scintilla- 
tion photons) below kinematic threshold [l3], having an 
effective cutoff at ~1 keV r . That both processes should 
decrease hand-in-hand can be argued based on the dom- 
inant role of ionization as the main precursor to direct 
scintillation for low energy recoils in LXe mil. 

• The energy scales postulated thus far by XENON10 
(color curves in Fig. 2) overestimate, by several orders of 
magnitude, the very small average charge yields observed 
in impact ionization experiments involving few keV and 
sub-keV xenon ions. The relevance of these measure- 
ments and examples from the literature are discussed in 
Q, where it is also emphasized that XENON 10 workers 
are not without a reference on what to expect at these 
low energies. The introduction of the adiabatic correc- 
tion proposed in [3(3 ] resolves this disagreement. 
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FIG. 2: Updated and adapted from ||. Color lines correspond 
to the rapid evolution of the XENON 10 S2 ionization energy 
scale over the past two years, using the "best-fit Monte Carlo" 
method [2j|. In chronological order: solid blue and green [2^ ]: 
dashed red [33l ] (workshop presentation); solid red [3oT ] (its 
proceedings); dotted red, la C.L. contour in [3^|, as per [{J; 
dashed blue [3, Efl] ■ Dark lines correspond to the predicted 
energy scales in [8| (see text). Inset: alternative analysis of 
the few-electron S2 background in XENON10 @] (see text). 



Of special interest is a population of XENON10 
low energy ionization signals described in [9j as single 
electrons 7 . This definition is both surprising and mis- 
leading. As mentioned in [8J, the large amplification gain 
provided by the electroluminescence provides a good res- 
olution in the multiplicity of drifted charge in LXe. These 
events clearly include a multi-electron component and 
their population has been described as such by the con- 
tact author of 9;], as recently as in (35[. Their origin is 
unknown, and hard to ascribe to minimum ionizing parti- 
cles in an efficient self-shielding medium such as LXe |8j . 
While none of this is discussed in Q , their accumulation 
towards low values of S2 pulse-width is to be expected 
from the effect of charge multiplicity on this variable, 
and does not have to correspond to a radioactive con- 
tamination close to the z = detector coordinate (this is 
unlikely, given that roughly the same number of PMTs, 
major sources of internal activity in the XENON10 de- 
tector, are placed at both extrema of z). 



7 The well-know ionization "afterpulses" following large energy de- 
positions in LXe, limited to single electrons but nevertheless 
mentioned in Q as a possible origin for these multi-electron 
events, can be trivially removed through a ~100/is delayed- 
coincidence cut 0, l2ll |34 | . No mention is made of the application 
of this cut to the dataset in [jj. This should be clarified. 
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IV: CONCLUSIONS 

The inset in Fig. 2 represents the differential rate of 
few-electron S2 events in XENON10, obtained by ap- 
plying the same five background cuts as in Q, extend- 
ing the analysis down to the S2 = le _ boundary, and 
adopting the energy scale described by the black solid 
line in the same figure. This differential rate seems to 
be also compatible with Fig. 2 in (35|, once the adopted 
energy scale is included. It offers a good match to the ex- 
pected signal from a light WIMP in the region of interest 
of other searches (DAMA/LIBRA, CoGeNT, CRESST). 
The outcome of this exercise, performed here strictly for 
the sake of argument, should be strongly de-emphasized 
at this time, given the present lack of knowledge about 
this energy scale evidenced in Fig. 2. The reader should 
remember instead that a further evolution by a mere ~1 
keV r in the ever-changing S2 energy scales postulated by 
XENON10 can transform the "severe constraints" of [9( 
into a signal in principle compatible with a light- WIMP. 

In conclusion, the claims in Q are clearly presently 
untenable. Awaiting clarification of the several pend- 
ing issues pointed out in Sec. I and II, light- WIMP lim- 
its obtained through a more conventional analysis of 
XENON100 data [lOj can only be assigned the very lim- 
ited meaning illustrated by Fig. 1. The XENON100 col- 
laboration is congratulated for the recent advancement 
in their understanding of £ e ff, encouraged to develop 
improved methods of characterization of the S2 energy 
scale leading to a reliable exploration of light- WIMP 
candidates, and urged to employ transparency in the 
discussion of uncertainties and assumptions underlaying 
their results, in view of the very limited performance of 
LXe as a light- WIMP detection medium. Finally, while 
several interesting phenomenological routes to alleviate 
tension between LXe constraints and other light- WIMP 
searches have been put forward recently [36| , these devi- 
ations from arguably more conventional assumptions do 
not seem to be mandatory at this time. 

N.B.: A new measurement of £ e ff by the ZEPLIN-III 
collaboration appeared coincident with the release of this 
preprint [37[. Fig. 1 now reflects XENON100 exclusions 
obtained with it. ZEPLIN-III derives a £ c g decreasing 
below 40 keV r and vanishing at few keV r , in tight agree- 
ment with the LXe kinematic threshold described in 1131 1 . 
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